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ITEM  IDENTITIES  AND  THEIR  RELATED  OBSERVABLES 


D.  P.  Gaver 
P.  A.  Jacobs 


Abstract 

Each  of  J  items  has  a  characteristic  Signature  which  varies  in  time.  At 
time  0,  the  value  of  a  Signature  and  the  identity  of  the  corresponding  item  are 
known.  No  further  values  of  Signatures  are  observed  until  a  later  time  t>0. 
At  time  t,  a  Signature  from  an  unknown  item  is  observed.  The  problem  is  to 
estimate  the  identity  of  the  item  whose  Signature  is  observed  at  time  t.  The 
estimation  procedure  studied  is  to  estimate  the  identity  of  the  unknown  item 
to  be  that  one  which  maximizes  the  posterior  probability  of  producing  the 
observed  Signature. 

Key  Words:  Classification;  Bayesian  Paradigm;  autoregressive  process 
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ITEM  IDENTITIES  AND  THEIR  RELATED  OBSERVABLES 


D.  P.  Gaver 
P.  A.  Jacobs 


1.  The  Problem 

We  are  concerned  with  a  diagnostic  problem  that  may  occur  in  many 
applied  areas,  with  medicine  and  mechanical  system  reliability  furnishing 
handy  examples. 

Suppose  there  are  J  Items  (e.g.,  diseases  or  faults),  each  of  which  has  a 
characteristic  Signature  which  varies  in  time;  the  Signature  of  Item  i  is 

Yi(t)  =  ei  +  Xi(t)  \Zoa%1.  (U) 

For  the  moment  {Xj(t)J  is  an  unspecified  univariate  stochastic  process,  but  one 
that  loosely,  has  the  property  of  staying  near  0j  in  finite  time,  and  for 
convenience  has  some  stationary  or  steady-state  behavior.  One  could  think  of 
Yj(t)  as  a  physical  index  characteristic  of  a  particular  disease  or  disability,  e.g., 
blood  pressure,  heart-beat  pattern,  mechanical  vibration  spectrum,  etc. 

Suppose  that  Yj(t)  is  only  evident  or  observable  occasionally  at  times 
unrelated  to  the  magnitude  of  Yj(t),  but  driven  by  other  forces  (this  is  a  bit 
implausible  in  some  situations  because  a  high  blood  pressure  or  temperature 
may  actually  induce  medical  observation;  the  formulation  can  be  changed  to 
reflect  this  later).  Suppose  that  the  Signature  and  the  identity  of  an  item  are 
both  observed  at  time  t=0.  Suppose  that  later  on,  at  time  t,  only  the  Signature 
of  an  item  is  observed.  The  first  question  is:  what  is  the  probability  that, 
given  the  Signature  value  observed,  its  originating  item  is  any  particular  one 
of  a  list  of  candidates,  e.g.,  i=13  out  of  a  particular  list  j=l,  2, ...  13, ...  J  =  39? 
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2.  Question:  Who  is  Being  Observed? 

Here  we  become  quite  specific,  assuming  for  illustration  that  the 
Signature  of  Item  j  is  AR(1): 


Yj(t)  =  0j  +  Xj(t), 


(2.1) 


where 


Xj(t)  =  pjXj(t-l)  +  Ej(t),  t  =  0, 1, 2, 17, ... 


(2.2) 


with  (£j(t)}  ~  HD  N  (0,Oj).  For  the  moment  assume  that  the  basic  constants  0j, 

2 

j  =  1,  2, ...,  J  are  known,  as  are  the  correlations  I  pj  I  <1  and  the  variances  Oj. 


If  I  Pjl  <1  then,  given  Xj(0),  Xj(t)  is  Normal/Gauss 


( 


Xj(t)  I  Xj(0)  ~  N 


XjW)pj,Oj(t)  =  Cj 


V 

or,  in  terms  of  the  actual  observable  Yj(t), 


1-p 


2t 


1-Pi 


\ 


J 


(2.3) 


Yj(t)  I  Yj(0)  ~  N((Yj(O)-0j)p|  +  0j,  o- (t)) 


Apparently  if  t  runs  on  indefinitely,  its  distribution  approaches 


(2.4) 


N^0j,Oj(oo)  j,  with  Oj(~)  =  oj/jl-pj  j 


(2.5) 


Scenario:  Suppose  there  are  potentially  J  items  around,  but  admit  the 
possibility  that  each  can  vanish,  so  that  there  is  some  probability,  which  may 
be  less  than  one,  that  the  item  will  actually  be  "around"  at  time  t.  Let  the 
probability  that  one,  the  jth,  is  present  at  t  be  pj  (t).  Next,  let  Sj(t)  be  the 
probability  that  the  jth  be  observable,  e.g.,  is  emitting  sounds,  given  that  it  is 
present;  assume  independence  between  items.  All  of  this  is  preliminary,  and 
can  (and  should)  be  modified  as  required.  Now  consider  that  event 

A)  Item  i  has  been  observed  at  t=0,  and  its  Signature  is  observed /noted  to 
be  Yi(0)  =  y(0). 
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If  it  has  been  a  long  time  since  an  observation  was  made  then  it  is 
plausible  to  take  the  probability  of  the  above  event  as 

P{ Yi(0)  £(dy(0)),  1(0)  =  i}  =  Pi(0)  Sj(0)(p  dy(0)  (2.6) 

(p  being  the  N(0,1)  density;  that  is,  we  observe  Item  i  after  his  Signature  is  in 
statistical  equilibrium,  and  measure  time  from  that  (such  an)  instant. 

Time  t  elapses  and  then  for  any  j  (either  the  same  as  i,  or  different) 

B)  Item  j  is  observed  at  time  t>0,  with  Signature  value  Yj(t)  =  y(t). 

Since  the  AR(1)  process  is  Markov,  we  can  write  down  the  joint 
probability 


P{Yj(0)  £(dy(0)),  1(0)  =  i,  Yj(t)  £(dy(t)),  I(t)  =  j}  = 


fy(0)-e^  dy(0)  _  fyW-mjft^dyft) 

Pi<°)  Si(0)  j-Jjfo-  Pi«>  ft)  ^  0|(l)‘  J  ^<0  (2-7) 


where  we  can  think  of  pj(t)  as  the  conditional  probability  that  j  is  actually 
eligible  for  observation,  and  Sj(t)  as  the  conditional  probability  of  being 
observable.  Note  that  it  is  natural  for  j*i  that 


mj(t)  =  &j 


(2.8,a) 


Oj(t)  =  Oj(oo)  = 


(2.8,b) 


for  at  t  it  is  still  a  "long  time"  since  j(*i)  was  last  observed;  whereas  if  j=i,  the 
same  item  that  was  observed  and  identified  at  t=0,  then 


and 


mj(t)  =  9*  +  (y(0) — 0j)  p| 


(2  9,a) 


(2.9,b) 


It  follows  that 
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C)  the  conditional  probability  of  observing  a  particular  item  and 
observing  it  with  Signature  value  y(t)  given  the  initial  condition  Cj(0)  = 
{Yi(0),I(0)}  is  just 


P{Yj(t)  e(dy(t)),  I(t)  =  j  I  Cj(0)  =  (y(0),  1(0)))  = 


m  m  ry(t)-mj(t)\dy(t) 

Pi(t)  si(t)<P[  Oj"( t)  J^(0 


(2.10) 


where  the  appropriate  mj(t)  and  Oj  (t)  come  from  (2.8)  or  (2.9).  In  words,  (2.10) 
provides  the  probability  of  seeing  or  positively  identifying  a  particular  item,  j, 
and  simultaneously  measuring  its  current  Signature  value,  given  such  a 
complete  observation  at  t=0. 


Next 


i 

D)  P{Y(0  £(dy(t))  I  Cj(0)}  =  ^  p{Yj(t)e(dy(t)),  I(t)  =  j  I  Q(0)} 


i 

=  IPj(t)Sj(t)  <p 


y(t)  -  mj(t)A 

.  <¥*>  J 


dy(t) 

tfj(t)‘ 


(2.11) 


This  is  the  probability  that  one  measures  a  Signature  at  t,  but  cannot  make  a 
complete  identification. 


From  (2.10)  and  (2.11)  we  can  calculate  the  probability  of  the  identity  of  the 
item  whose  signature  is  measured  at  time  t,  given  the  last  complete 
observation: 


E)  P(I(t)  =  j  I  Y(t)  =  y(t),  Cj(0)}  = 


Pj(t)  Sj(t)  (pj 

ry(t)-mi(t)> 

,  °i(1)  J 

1 

Oj(t) 

1 

Jpk(t)  Sk(t)  <P 

^y(t)  —  mk(t)> 

1  tfk(t)  J 

1 

Ok(t) 

(2.12) 


If  one  looks  hard  at  (2.12)  it  seems  clear  that  when  pj>0  and  y(t)  is  near 
y(0),  and  j=i,  the  above  probability  becomes  closer  and  closer  to  unity  as  t 
decreases,  while  if  j*i  that  probability  (that  Item  j#i)  correspond¬ 
ingly  decreases. 
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3.  The  Probability  of  Making  the  Wrong  Decision 

In  this  section  we  study  properties  of  making  a  decision  based  on 
choosing  the  item  which  maximizes  the  posterior  probability  of  a  particular 
item  being  the  one  whose  Signature  is  observed  at  time  t  given  the  last 
complete  observation  at  time  0.  For  simplicity  we  will  assume  that  there  are 
two  items  whose  parameters  are  0i  and  02- 

3.1  The  parameters  9j  and  6^  are  known. 

Suppose  the  parameters  0]  and  02  are  known  and  I  Pil  <1.  The  condi¬ 
tional  distribution  of  Y(t)  given  Y(0)  =  y(0),  1(0)  =  1,  I(t)  =  1  is  normal  with 
mean 


mi(t)  =  0i  +  (y(OMh)p* 


and  standard  deviation 


oi(t)  =  ojM 


Further, 

Y(t)-m2(t)  Y(t)-mi(t)  +  (mi(t)-m2(t)) 
C2(t)  =  02(0 

Ql(t)  Y(t)-mi(t)  mi(t)-m2(t) 
"  02(t)  Gl(t)  +  02(0 


(3.1) 


(3.2) 


(3.3) 


where  n\2(t)  and  mi(t)  are  given  by  (2.8). 

Assume  that  item  j  will  be  identified  with  the  Signature  at  time  t  if  it 
maximizes  (2.12);  then. 


P  {wrong  item  is  identified  1 1(0)  =  1, 1(t)  =  1,  Y(0)  =  y(0)}  (3.4) 


1  (Y(t)-m2(t))2 

► 

=  F" 

p2<t)s2(t)^-02(t)  exp 

“2  2 
a2(t) 

r  ; 

j] 

1 

1  (Y(t)-mi(t)) 

- 

>p,(,)s,(t)  Vtoo,(o  exp 

2  oi(t)2 

> 

I(0)=1, 1(t)=l,  Y(0)=y(0)  f 


j 
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■  =■  exp  f! z2 + \  (cMZ + d<t))2] 


where  Z  has  a  standard  normal  distribution 


c(t)  "  c2(t)  " 


^  1-Pi‘ 


ciO”) 


d(t)  = 


mi(t)  -  m2(t)  (y(0Hl)t)>  +  (01 "  02) 


O2(oo) 


P  (wrong  decision  |  1(0)  =  1, 1(t)=l,  Y(0)  =  y(0)} 


2(t)s2(t) 


,  2t  ,  . 
l~Pl  01  (oo) 


pi(t)si(t)  o2(«0 


expj-^  [(l-^c(t)2)Z2  -  2c(t)d(t)Z  -  d(t)2  ] 


=  P  -2  In 


2(0s2(t) 


f  2t 

1-Pl  Ol(°°) 


pi(t)si(t)  02(°o) 


<  (l-c(t)2)Z2  -  2c(t)d(t)Z  -  d(t)2 


=  p{o<(W(.)2)z2-2c«»d(t)Z-d(.>2  4  2ln  } 

Note  that  as  t->0  lnj^l-p^j-^-oo  and  thus 


and  thus 


P{wrong  decision  I  1(0)  =  1, 1(t)  =  1,  Y(0)  =  y}  — ►  0. 

Example:  Suppose  Oj  =  a  i  =  1,2,  1  pi  I  =  I  p  I  <1,  i=l,2  and  pi(t)si(t)  =  p2(t)s2(t); 
then 

p{wrong  decision  I  1(0)  =  1, 1(t)  =  1,  Y(0)  =  y(0)} 

=  P  {0  <  (1— c(t)2)Z2  -  2c(t)d(t)Z  -  d(t)2  +  2  In  c(t) }  (3.7) 
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with 


c(t)  =  V 1— p2t 


(3.8,3) 


and 


,,  „  (y(0)  -  ©Op*  +  (61-62) 

d<t,=  “aTVT^ 


(3.8,b) 


(01-02) 

As  t— d(t)-» —  7==-7  c(t) — >1  and 


o/Vl-p2 

P{wrong  decision  I  1(0)  =  1, 1(t)  =  1,  Y(0)  =  y(0)} 

i  Id(~>'  " 

1-0 


(3.9) 


where  O  is  the  standard  normal  distribution  function.  In  this  case,  the 
conditional  probability  of  making  a  wrong  decision  decreases  as  I  0j— 02  I 
increases  and  the  probability  increases  as  a  increases. 

For  arbitrary  t  >  0  such  that  c*l, 

P{wrong  decision  I  1(0)  =  1, 1(t)  =  1,  Y(0)  =  y(0)} 


=  Pl 


r  c(t)  d(t)> 

2  d(t)2  2  In  c(t)| 

~~  1— c(t)2  , 

1  >(l^(t)2)2  1_c(t)2] 

_ i 

„  c(t)  d(t) 

/  d(t)2  2  In  c(t) 

Z  1— c(t)2 

\  (l^(t)2)2  1_c( 

„  c(t)  d(t) 

/d(t)2  In(l-p2t) 

► 

>  \  p4*  p2t 

=  1  -  <D  (  Vl-p21  p~2td(t)  +  k(t))  +  O  (  Vl-p2t  p~2td(t)  -  k(t))  (3.10a) 


where 


k(t)  =  [d(t)2  p-4*  -  p~2t  In  (l-p2t)j  2 . 


Note  that  as  t— >0,  c(t)->0  and  the  conditional  probability  of  an  incorrect 
decision  tends  to  zero. 
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A  numerically  more  stable  form  of  (3.10a)  is  the  following. 


P{wrong  decision  1 1(0)=!,  I(t)=l,Y(0)=y(0))  =  l-4>(Ai)+<t>(A2/  (3.10b) 


where 

A!  =  [-d(t)2  +  21n  c(t)]  [c(t)d(t)  -  y/d(t)2-  (l-c(t)2)21n  c(t)]'1 

A2  =  [— <d(t)2  +  21  n  c(t)]  [c(t)d(t)  +  Vd(t)2- (l-c(t)2)21n  c(t)] 

The  unconditional  probability  of  a  wrong  classification  can  be  found  by 
taking  the  expected  value  of  (3.10b)  with  respect  to  the  normal  random 
variable  Y(0)  with  mean  «i  and  standard  deviation  a(°°).  The  evaluation  of 
the  integrals  is  performed  numerically  using  Gauss-Hermite  quadration  with 
59  points;  cf.  A.  H.  Stroud  et  al  [19661  and  Naylor  and  Smith  [1982],  Figures  1 
and  2  show  the  unconditional  probability  of  wrong  classification  at  time  1,  2, 
...,  100  for  parameter  values  0i  =  1,  82  =  2,  a  =  1  and  values  of  p  =  0,  0.5,  0.8,  0.90, 
0.95.  (resp.  curves  A,  B,  C,  D,  E). 

Figure  1  shows  the  probability  of  wrong  classification.  As  noted  above, 
the  probability  of  wrong  classification  is  small  for  small  values  of  t  and  then 
increases  to  a  limiting  value.  For  very  small  times  having  a  larger  p  results 
in  smaller  probabilities  of  error.  It  is  seen  from  (3.9)  that  the  limiting  value  is 

a  function  of  o(°°)  =  o/V  1-p2  .  Hence  the  larger  p  is  the  larger  the  limiting 
probability  of  wrong  classification.  Thus,  while  for  very  small  times,  having  a 
large  positive  p  decreases  the  probability  of  wrong  classification,  at  larger  time 
it  increases  the  probability  of  wrong  classification. 

In  Figure  2,  each  probability  of  wrong  classification  is  divided  by  its 
asymptotic  value.  These  normalized  probabilities  are  between  0  and  1.  The 
fraction  (relative  to  the  limiting  probability)  of  the  probability  of  wrong 
classification  for  each  time  t  is  monotone  in  p.  For  a  fixed  time  t,  for  a  high  p, 
the  fraction  of  the  probability  of  wrong  classification  is  less  than  or  equal  to 
that  for  a  lower  p.  The  higher  p  is,  the  longer  it  takes  for  the  probability  of 
error  to  reach  its  asymptotic  value. 

3.2.  The  parameters  0i  and  62  have  noninformative  priors. 

In  this  subsection  we  will  assume  that  0]  and  02  have  noninformative 
prior  distributions.  Thus,  assuming  I  Pi  I  <1, 
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PROBABILITY  OF  WRONG  CLASSIFICATION 
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Figure 


NORMALIZED  PROB(WRONG  CLASSIFICATION) 
TH1  =  1  ;TH2=2;SIG=1 
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Figure 


(3.11) 


p{y(0)  edy(O),  eedel  i(o>  =  1} 

=  -== -  exp  -7  (y(0)  -  0)2/oi(««)2  dy(0)d6 . 

\  2n  01  (“°)  z 

Hence,  formally  for  a  non-informative  improper  prior, 

P  {Y(0)edy(0)  1 1(0)  =  l)  =  l-dy(O) 
and  thus 

P  {eed0l  1(0)  =  1,  Y(0)  =  y(0)} 

=  p-— exp  \-l  (e-y(0))2/oi(~)2  d0  . 

\  27101  (°°)  1 

Further, 

P  (Y(t)edy(t),  0£de  1 1(0)=1, 1(0=1,  Y(0)=y(0)} 

1 

V2 n  oi(t) 

1  1  2  2 

x  -7= — 7 —  exp  -p  (6  -  y(0))  /o iW 
y2n  <J](°°)  ^ 

1  1  2  2  1 

=  -== - — -  exp  -5 (0-pi(t))  /vi(t)-~R 

V27toi(t)V2^ai(oo)  >  2 

where  the  usual  calculations  yield 


exp  |-|  ^y(t)  -  0-p  1  (y(0)  -  9)  J /oi(t) 


(3.12) 


(3.13) 


(3.14) 


(3.15) 


(3.16) 
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R  = 


vi(t) 

ai(~)  oi(t) 

.2  2 


[y(t)  -  y(0)]2  =  [yCO-yCO)]2 /2 


l-Pi 


Cl(oo) 


Ol(oo)  =  O?  /  ^1-p  j  j 

Ol(t)2=  ^1-pf  joi(oo)2 
Hence, 

P  {v(t)edy(t)  I  1(0)  =  1, 1(t)  =  1,  Y(0)  =  y(0)} 

/ 

1 


1-p‘l 


exp  |-|(y(t)-y(0))2/  2  1— p\  ai(°°) 2 


Similarly 

p{Y(t)edy(t),0ed0 1 1(t)=2, 1(0)=1,  Y(0)=y(0)} 
=”7"7  exp|-l  (y(t)-0)2/o2(°°)2  dy(t)d0 

Ooloo)  l  *• 


V  27C  02(  oo) 


and 


P  (Y(t)edy(t)  I  I(t)=2, 1(0)=1,  Y(0)=y(0)}  =  l-dy(t) 
where  02(°°)  =  02/  fl-p2j  • 

Suppose  item  j  is  estimated  to  be  the  one  observed  at  time 
maximizes 


I(t)=j 


1(0) 


=1,  Y(0)  =  y(0),  Y(t)=y(t)}  . 


It  follows  that 


P  {item  2  is  chosen  I  I(0)=1,  Y(0)=y(0),  I(t)=l} 


(3.17) 

(3.18) 

(3.19) 

(3.20,a) 


(3.20,b) 

t  if  it 
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=  P  Ip2(t)s2(t)  >  pi(t)si(t)  ~f=i — 7-  expf-5  (V(t)— y(0))2/<X](t)2  }  1 I(0)=1, I(t)=l,  Y(0)=y(0)] 

V  2n  ai(t)  l  1  I 


=  P{^T aiW  yj  2jVp^  >  exp{-|  Z 

V~2n 


=  P. 


-2  In 


2 

Pi 


pi(t)si(t) 


<Z 


(3.21) 


where  Z  is  a  standard  normal  random  variable  and  ai(t)  =  fl-pjj  Oi(°°). 


Note  that  as  t— >0,  the  probability  of  choosing  the  wrong  object  tends  to  0. 
If 


k(t)  >  VaT  a/ 2(l-p‘l)  >1  (3.22) 

\1-Pl 

then  the  conditional  probability  of  making  the  wrong  decision  is  1. 

If  k(t)  <  1  then 

P  {item  2  is  chosen  |  I(0)=1,  Y(0)=y(0),  I(t)=l} 

=  2  [l-<D(-2  ln(k(t)»]  (3.23) 

where  <I>  is  the  standard  normal  cdf.  Note  that  (3.23)  is  independent  of  y(0). 
As  ai  increases,  the  conditional  probability  of  making  the  wrong  decision 
increases.  Similarly  as  t  increases  the  probability  of  making  the  wrong 
decision  increases.  As  pi(t)si(t)  decreases,  the  probability  of  making  the  wrong 
decision  also  increases. 


In  the  case  in  which  pj(t)  si(t)  =  p2(t)  S2(t),  pi  =  P2  =  P>0,  and  Oi  =  02  =  o 
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For  t  =  1,  k(t)  and  hence  the  probability  of  error  decreases  as  p  increases.  For 

t  =  2,  k(t)  =  a  for  all  values  of  p  and  so  the  probability  of  choosing  the 
wrong  item  increases  as  a  function  of  t.  The  limiting  value  of  k(t)  as  t-»«> 

once  again  depends  on  c(°°)  =  al^j  1-p2  .  Hence  as  p  increases  the  limiting 
probability  as  t-»«>  of  choosing  the  wrong  item  increases.  Figure  3  displays  the 
probability  of  choosing  the  wrong  item  for  a  =  0.1  and  p  =  0,  0.5,  0.8,  0.9,  0.95 
(resp.  curves  A,  B,  C,  D,  E).  Note  that  for  t>2,  it  is  p  =  0  that  has  the  smallest 
probability  of  incorrect  identification.  These  results  depend  on  the  specific 
prior  (3.20,b).  In  the  next  subsection  we  will  consider  Gaussian  priors. 

3.3  The  parameters  0i  and  62  have  Gaussian  priors 

In  this  section  we  will  assume  0i  has  a  normal  prior  with  mean  pi  and 
variance  ij  ,i  =  1,2.  Thus,  assuming  I  pi  I  <1, 


P  {Y(0)edy(0),  0ed0 1  l(0)=l} 


V  2tc  01(00)  V  2k  x\ 


exp 


1  (y(0) — 0)2  1  (8-pi)21 

'2"o^T'2  t2 


dy(O)d0 


(3.24) 


V  2n  ai(  00  )  V  2jc  T] 


exP| 


1  (0-m,(O))  1 

■5  1 —  -9  Ri(0) 

2  V,(0)  2 


dy(O)d0 


where  straightforward  calculations  yield 

,  ^2  2 

-2  1  1 

vi(0)  =  ~2  +  —J  =  "T 


T 


X\  oi(«)  ti  Cl(oo) 


mi(0)  =  (y  (0)  /  01  (<x>)2  )  +  jVi/xfj 


vi(0) 


2  2 
X\  y(0)  +  01  (00)  m 

;  2 

01  (00) 


(3.25) 


(3.26) 


(3.27) 
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Rl(0)  = - 2 - T  (y(0)  -  (ii)2 

ai(oo)  +  xi 


with 


(3.28) 


ai(oo) 


Therefore, 


P{eed0|l(O)=l,Y(O)=y(O)}  = 


V  2n  vi(0)  exp 

Similarly,  from  (2.4) 

p{Y(t)edy(t),  0ed6 1 1(0)=1, 1(t)=l,  Y(0)=y(0)} 


1  (0-mi(O))^ 


V^2tc  oi(t) 


exp 


(y(t)  -  0  -  Pi(y(O}-0))2/ai(t)2 


1  1  2  2) 

X  ^J~2k~v\(0)  eXpi“2  / vi(0)  dy(t)d0 


V  271  ai(t)  V  2ti  vi(0) 


exp 


1  (&-ai(t))  i 
"2  ^“^2 —  — 2  Ri(0  j"  dy(t)d0 


where 


^i(t)  =  (l-pi)2/oi(t)2  +  vi(0)  2j 


(l-p',)[y(t)-p*y(0)]  mi(0) 


«l(t)  =  ^l(t) 


+  ~2 
ai(t)  vi  (0)  J 


(3.29) 


(3.31) 


(3.32) 


(3.33) 


and 


(3.30 
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r  -p  p  (t)2 

Rl(t)  =  (y(t)-pt1y(0))-mi(0)(l-pt1)  — *5 - 3- 

J  ai(t)  vi  (0) 

=  (y(t)-pt1y(0))-mi(0)(l-pt1)  2/  (l-pt1)2vi(0)2+  ai(t)2  j 


2  2  2t  2 

oi(t)  =  Oj  (1— pi  )/(l— pi). 


Thus, 

P  (Y(t)edy(t)  I  Y(0)=y(0),  I(0)=1, 1(t)=l} 

‘  V2nt,(.)  6Xp{^(y<,)-W<,))2/T,(t>2 


where 


ti(t)2  =  (l-pi)2vi(0)2  +  oi(t)2 


=  ai(«) 


(l-p2t)ai(oo)2+  x2  2(1-^) 


ai(oo) 


=  (1— Pi)  ai(«>) 


1  +  T2  +  p\  O\(o»)2  /jx2  +  ai(oo)2  j 


(3.34) 


(3.35) 


(3.36) 


Hl(t)  =  piy(0)  +  mi(0)  (1-Pi) 


(3.37) 
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p}y  (0)1  o i  (°°)  +X\  +  Xiy  (0)  +  oi  (<*)  m  (1  -  pi) 


ai  («>)  +  x  i 


x]  y(0)+oiM2jii  j+pi  y(0)|ai(«>)2  +  Tj  j  -  x  jy(0)  +  Oi(oo)2p1 


Cl(oo)  +X} 


2  2  t  2 

Xj  y(0)  +  ai(«>)  Hi  +  P  iOj(°°)  [y(0)-pi] 

;  j 

ai(~)  +t! 


Hioi(~)2(l-p  i)  +  y(0)  x2  +  plai(oo)2 


Ol(oo)  +  Xl 


Similarly, 


P  (Y(t)edy(t)  I  Y(0)=y(0),  I(0)=1, 1(t)=2}  =  exp 

V  2ji  x2(t) 


1  (y(t)-H2(t>y 

2  5 


where 


2  2  2 
X2(t)  =  02(0°)  +  X  2 


H2(t)  =  H2- 


(3.39) 


(3.40) 


Once  again  we  will  estimate  the  identity  of  the  item  that  was  observed  at 
time  t  to  be  that  item  j  for  which 


(t)Sj(t)  1,  ,  ^2,  /s2) 

Ti(t)  "  exP  -2  (y(t)-H(t))  /xi(t) 


(3.41) 


is  the  largest. 

Thus, 

P  {choose  item  2 1  Y(0)=y(0)  ,  I(0)=1,  I(t)=l} 


=  P  {  &X2(0  *  eXp  { ^Y(t)  “  ^2(t))2/T2(t)2  J 

pl(t)si(t)  1  2  2  1, 

>  r-Tl(t)  exp  -  2  (Y(tHii(t))  /ti(t)  I  Y(0)=y(0),  I(0)=1, 1(t)=l 

[  P2(t)s2(t)  j  1  w_2)  pi(t)si(t)  f  1„21  1 

=pi  expi  -2(c<t)z+d(t»  J>  exp i ~2Z  } 


(3.42) 


(3.43) 


where  Z  is  a  standard  normal  random  variable. 


c(t)  = ; 


(3.44) 


oiW'y  i-pi 
a /a 2M  +t2 


*1  +  piaiM2Vfxi+  oi(oo)2!  1 2 


d(t)  = 


M-l  (t)-M-2(t) 


'  7^2  2 

02\°°)  +  ^2 


Hence 


(3.45) 


t  2  2  2 

p](y(0)-pi)  ai(oo)  +  (y(0)  -  P2)  1 1  +  (P1-P2)  ai(«») 


CTl(oo)  +Tj 


P  {choose  item  2 1  Y(0)=y(0)  ,  I(0)=1,  I(t)=l} 


(3.46) 
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d(t)2 

l-c(t)2 


-2  In 


p2(t)S2(t) 

pi(t)si(t)_ 


< 


Vl-c(t)2Z- 

V 


c(t)d(t)  > 

2 

r-  2 

- 

\l-dt)  J 

if  I  c(t)  I  <1 


If  pi(t)  si(t)  =  0,  then  item  2  will  always  be  chosen.  Assume  pi(t)  silt)  >  0.  < 

If  p2(t)s2(t)=0/  then  clearly  item  2  will  never  be  chosen.  In  what  follows 
assume 

pi(t)  si(t)  and  p2(t)  S2(t)  are  both  positive.  Note  that  as  t-»0,  Ti (t) — >0,  c(t) — >0 
and 


P  {choose  item  2  I  Y(0)=y(0)  ,  I(0)=1,  I(t)=l}  ->  0. 

Also,  as  and/or  O2(c»)->oo/  c(t) — > 0,  d(t) — > 0  and  the  conditional 

probability  of  making  a  wrong  decision  tends  to  zero  for  t>0.  Note  that  as 
ai(t)— c(t) — d(t)-»«>,  and  the  conditional  probability  of  making  a  wrong 
decision  tends  to  1. 

The  unconditional  probability  of  a  wrong  classification  can  be  found  by 
taking  the  expected  value  of  (3.46)  with  respect  to  Y(0).  It  follows  from  (3.25) 
that  Y(0)  has  a  normal  distribution  with  mean  pi  and  variance  Gj  (°°)2  +  X}2. 
Figures  4  and  5  show  the  probability  of  misclassification  for  the  case 
P2(t)  S2(t)  =  Pi(t)  sj(t),  ti  =  x2  =  1,  G]  =  02  =  1,  Pi  =  1,  P2  =  2  and  p  =  0.0,  0.5,  0.8, 
0.9,  0.95  (resp.  curves  A,  B,  C,  D,  E).  The  integration  to  evaluate  the  probability 
is  performed  using  Gauss-Hermite  quadration  with  59  points;  ef.  A.  H.  Stroud 
et  al.  [1066]  and  Naylor  and  Smith  [1982].  Figure  4  shows  the  probability  of 
wrong  classification.  For  small  times  t  having  a  larger  p  results  in  smaller 
probabilities  of  error.  The  probability  of  wrong  classification  is  small  for  small 
values  of  t  and  then  increases  to  a  limiting  value  as  t->°°.  The  limiting 
probability  of  wrong  classification  increases  as  p  increases.  Thus,  while  for 
small  times,  having  a  larger  positive  p  decreases  the  probability  of  wrong 
classification,  at  larger  times  it  increases  the  probability  of  wrong 
classification.  In  Figure  5,  each  probability  of  wrong  classification  is  divided 
by  its  asymptotic  (as  t— >°°)  value.  For  a  fixed  time  t,  for  a  higher  p,  the  fraction 
of  the  probability  of  wrong  classification  is  less  than  or  equal  to  that  for  a 
lower  p.  The  higher  p  is,  the  longer  it  takes  for  the  probability  of  error  to 
reach  its  asymptotic  value. 
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Fi  gure 


4.  Conclusions 

This  paper  considers  a  model  in  which  there  are  J  items  each  of  which  has 
a  characteristic  Signature  which  varies  in  time.  In  the  examples  studied  the 
Signature  process  of  each  item  is  an  AR(l)-type  process  independent  of  the 
other  items.  At  time  0  the  value  of  a  Signature  and  the  identity  of  the 
corresponding  item  are  known.  No  further  values  of  Signatures  are  observed 
until  a  later  time  t.  At  time  t  >  0,  a  Signature  from  an  unknown  item  is 
measured.  The  problem  is  to  estimate  the  identity  of  the  item  observed  at 
time  t.  The  estimation  procedure  studied  is  to  estimate  the  identity  of  the 
unknown  item  to  be  that  one  which  maximizes  the  posterior  probability  of 
producing  the  observed  signature.  The  probability  of  correct  identification 
depends  on  the  parameters  of  the  processes  and  the  magnitude  of  t.  The 
smaller  t  is  the  greater  the  probability  of  correct  identification  of  the  unknown 
item.  For  fixed  time  t  a  higher  positive  correlation  results  in  a  lower  fraction 
of  the  probability  of  wrong  classification  with  respect  to  its  limiting  (as  t— >«») 
value. 

It  appears  that  the  problem  of  this  paper  resembles  the  "missing-data 
problem"  of  statistics;  see  Little  and  Rubin  (1987).  We  plan  to  investigate  the 
connection  more  extensively  in  future  work. 
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